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NOTATION 


a 

Y 

C D 

C L 

c e 

D 

d i 

L 

*i 

R 

s i 

T 

V 

W 

X 


x 

y 

z 


z 

Az 


total acceleration, m/sec 2 
longitudinal acceleration, m/sec 2 
normal acceleration, m/sec 2 
lateral acceleration, m/sec 2 
drag coefficient 
lift coefficient 

d i “ s i ( see fig- 2 ) 
aerodynamic drag, N 

ground- track distance between way points 
aerodynamic lift, N 

direction cosine of ith ground- track segment, 
ith ground- track segment length, m 
turn radius in ground track, m 
ground- track distance from last way point to next turn (see fig. 2) 
thrust, N 

longitudinal velocity, m/sec 
weight of vehicle, N 

/ x i\ 

position vector in horizontal plane,/ J, m 

\Yi/ 

displacement in x direction of local earth-fixed coordinate system, m 
displacement in y direction of local earth-fixed coordinate system, m 

/ a i\ 

position vector in vertical plane, j j, m 

displacement in z direction of local earth-fixed coordinate system, m 
change in z as a result of vertical acceleration (see fig. 3) 
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a angle of attack, deg 

Y flight path angle, angle between velocity vector and horizontal plane, deg 

A distance along the trajectory, m 

6j? flap deflection, deg 

n thrust angle, between body axis and thrust vector, deg 

0 change in inclination angle, deg (see fig, 4) 

0^ change in heading of jth curved segment, deg 

p radius of flight path in vertical plane, m 

a distance along ground track, m 

4> roll angle, angle about the velocity vector, deg 

heading angle, angle between a reference direction in horizontal plane 
and projections of V in horizontal plane, deg 

Note: Italics for the symbol denote quantities defined in a vertical plane 

and are analogous to those nonitalicized letters in the ground-track plane. 
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OPTIMAL GUIDANCE AND CONTROL FOR INVESTIGATING 


AIRCRAFT NOISE- IMPACT REDUCTION 
Elwood C- Stewart and Thomas M. Carson 
Ames Research Center 


SUMMARY 


As part of a NASA program to investigate technical approaches for increas- 
ing the terminal area effectiveness of advanced short-haul aircraft, this 
report is concerned with a methodology for investigating the reduction of com- 
munity noise impact. There are a number of computer programs available for 
the generation of noise footprints, all of which require various input data 
that are usually unavailable. This report is concerned with the development 
of two models to provide such data: a guidance generator and an aircraft con- 

trol generator suitable for various current and advanced types of aircraft. 

The guidance generator produces the commanded path information from inputs 
chosen by an operator from a graphic scope display of a land-use map of the 
terminal area. The guidance generator also produces smoothing at the junc- 
tions of straight-line paths. The aircraft control generator determines the 
optimal set of the available controls such that the aircraft will follow the 
commanded path. The solutions for the control functions are given and shown 
to be dependent on the class of aircraft to be considered, that is, whether 
the thrust vector is rotatable and whether the thrust vector affects the 
aerodynamic forces. For the class of aircraft possessing a rotatable thrust 
vector, the solution is redundant; this redundancy is removed by the addi- 
tional condition that the noise impact be minimized. Information from both 
the guidance generator and the aircraft control generator is used by the foot- 
print program to construct the noise footprint. 

The complete package of models provides a useful methodology for con- 
structing an interactive graphics tool for rapidly assessing the effects of 
different aircraft technologies, flight paths, aircraft mixes, and other vari- 
ables, and for the minimization of noise impact. These effects are illus- 
trated, utilizing an existing footprint program, by several examples which 
consist of three broad classes of CTOL and V/STOL aircraft. 


INTRODUCTION 


NASA has been engaged in a program to investigate technical approaches 
for increasing the terminal area effectiveness of advanced short-haul aircraft. 
The program is intended to provide guidance to NASA’s advanced aircraft tech- 
nology program. The principal purpose of the program is to investigate the 
impact of various technical alternatives for increasing the effectiveness of 
short-haul aircraft in terms of airport airside capacity and delay, community 



noise , air pollution, fuel consumption, and ride quality. The technical 
alternatives include various advanced aircraft, their greater operational 
bounds and range of flight paths, and avionics. 

This report is concerned with the effects of aircraft noise. Because of 
public concern with environmental issues, the importance of the effects of 
aircraft noise has increased significantly in recent years. We are concerned 
here with developing a methodology for examining the trade-offs between tech- 
nical alternatives and community noise impact. 

A key tool in the assessment of noise impact has been the development of 
noise footprint models and programs. Common to all these models and programs, 
however, is the need to generate reliable input information so as to render 
the results more realistic and to enhance their use in examining the many 
variables and alternatives associated with noise reduction. 

The input information is of two distinct types. First, a noise footprint 
program must be supplied with a description of the aircraft path. For three- 
dimensional path optimization purposes, the desired aircraft path is best 
supplied by a guidance generator which utilizes basic information from land- 
use maps. In this way, the guidance generator generates desired or optimal 
paths which are directed over areas that are less sensitive to noise. The 
guidance generator also generates smoothing, that is, smooth commanded paths 
at the junctions of straight-line segments so as to avoid unrealistically high 
accelerations . 

Second, the noise- footprint program must be supplied with aircraft con- 
trol function information in order to determine the noise footprint based on 
stored information in the footprint program. How much of the control informa- 
tion is required depends on the particular footprint program. Most programs 
inherently assume isotropic noise characteristics, and for them only one of 
the aircraft controls, the engine parameter (such as engine thrust, throttle 
setting, etc.) is required. A more sophisticated footprint program that could 
account for anisotropic noise characteristics would require all the aircraft 
control functions — the engine parameter, the directional control of the 
active thrusting forces, bank angle, and angle of attack. In any case, it is 
necessary to relate the aircraft control functions to the path being flown as 
specified by the guidance generator. This is best accomplished by utilizing 
the aircraft equations of motion to derive an aircraft control generator that 
determines the optimal set of the available controls from which those required 
by the footprint program are taken. 

In this report, we derive the models for generating the input information 
required by noise footprint models or programs. The two models are a guidance 
generator and an aircraft control generator suitable for various types of cur- 
rent and advanced aircraft. These two models will be made compatible with a 
particular noise footprint program, widely used in Government and industry, 
that was developed several years ago by Serendipity for the Department of 
Transportation (refs. 1 and 2). The models could be used with a more general 
footprint program taking into account anisotropic noise since all the required 
information is available from the models developed here. 
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The complete package, when implemented, provides a useful interactive 
graphics optimization tool for rapidly assessing the effects of different air- 
craft technologies, flight paths, aircraft mixes, and other variables, and for 
minimizing of noise effects. Several examples are given to illustrate appli- 
cations of the methodology to broad categories of advanced aircraft. 


GENERAL STRUCTURE 


The general structure of a system incorporating a noise footprint model, 
the guidance generator, and the aircraft- control system model are shown in 
sketch (a) . 


STORED 

INFORMATION 



ACTUAL TRAJECTORY 


Sketch (a) 

The guidance generator produces the desired trajectory and the control system 
then attempts to make the aircraft follow that trajectory. The final element 
is the footprint program which requires three kinds of information: 

(1) stored information describing the noise characteristics of the aircraft; 
for the footprint program used here this stored information consists of the 
noise characteristics expressed as a function of the slant range and the 
engine noise parameter; (2) input information describing the aircraft trajec- 
tory; and (3) input information defining the variation of the aircraft con- 
trols during the flight; for the footprint program to be used here, only the 
engine parameter is utilized. 

It will be desirable to modify the structure from that just discussed to 
accommodate the requirements of the particular noise footprint program to be 
used, and to achieve certain simplifications. In this regard, there are three 
considerations . 

First, the footprint program requires the trajectory to be parameterized 
in terms of segment lengths, parts of circles, and inclination angles. 

Second, since the aircraft can follow the commanded trajectory exactly (in the 
absence of control errors), the system structure shown above can be simplified 
as shown in sketch (b) on the following page. 

Third, as a result of the modified structure, it is necessary to generate 
the optimal aircraft controls based on the desired trajectory, as indicated in 
the sketch by the aircraft control generator. The situation shown is for the 
footprint program to be used here in which only one of the controls, the 
engine parameter T is used. For a more general program, all the controls 
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Sketch (b) 

are available for use. As will be seen later, the additional controls consist 
of the directional control of the active thrusting forces n» bank angle cf), 
and angle of attack a. This information, in conjunction with the trajectory 
information, provides all the information that would be required by any foot- 
print program. 


GUIDANCE GENERATOR 


Introduction 

The purpose of the guidance generator is to generate information about 
the desired path of the aircraft. In this section we will consider how to 
generate this information from available data and how to generate it in a form 
that will be compatible with the existing footprint program. 

First, consider the input information available to the guidance gener- 
ator. The determination of the desired or optimal ground track information 
for the purpose of minimizing noise effects must necessarily start with land- 
use maps of the airport and the surrounding community. This will permit the 
selection of candidate flight paths that are directed away from noise- 
sensitive areas. The most convenient way in which to display such maps is on 
a graphic scope display since this achieves good man-machine interaction for 
the noise minimization process. In this form, the graphics display readily 
incorporates provision for the experimenter to specify certain trial way- 
points that are the start and end points of straight-line segments of the 
ground track. The coordinates of these way-points are readily stored by the 
computer and hence are available as inputs to the guidance generator. To 
completely define the path of the aircraft it is also necessary to specify the 
velocity and inclination angle along each segment. And finally, it will be 
necessary to specify the tolerable acceleration levels to enable the genera- 
tion of suitably smoothed transitions between straight-line segments, as will 
be discussed shortly. 

Next, consider the form of the output of the guidance generator. Since 
the guidance generator output will be used directly by the existing footprint 
program which is to be left unchanged, this output must be of a form suitable 
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for the footprint program. The latter requires that the flight path be 
described by a sequence of segments in which each segment is defined by the 
following four variables: 

or 0^. kth segment length in horizontal plane, m, or kth subtended turn 
angle, deg 

y inclination angle, climb or descent, for kth segment, deg 

K 

engine noise parameter for the kth segment, percent of maximum 

R^. turn radius of kth segment, m 

0 = straight 
+ = clockwise 
- = counterclockwise 

The parameters, 1 %, 0^, and R^, need to be provided by the guidance generator, 
the Yfc is specified as an input, and the T^ is to be determined by the 
aircraft control generator. 

An important consideration in the generation of the commanded trajectory 
is the need to generate a smooth transition between straight-line segments so 
as to avoid unrealistic situations at the junction of straight-line segments. 
This will be done in the general three-dimensional case. 

The transition between straight-line segments is characterized by various 
combinations of lateral, normal, and longitudinal accelerations. There are 
seven possible combinations of motion on the transition segments as enumerated 
by the following: 


Case 

• 

V 

* 

Y 


1 

0 

^0 

^0 

2 

0 

0 

^0 

3 

0 

#0 

0 

4 

^0 

0 

0 

5 

^0 

0 

^0 

6 

^0 

^0 

0 

7 

^0 

*0 

#0 


For each of these cases, it will be necessary to generate the information 
required by the footprint program. However, many of the cases are similar and 
the results for some may be constructed by combining other cases, as will be 
seen in later developments. Basically, there are three types of transition 
trajectories which need to be considered: 

Type A: arcs with free end points 

Type B: arcs with constrained end points 

Type C: longitudinal acceleration with fixed end points 

Rather than discuss each of these types in general and then apply them to each 
of the seven possible cases, it will be more meaningful to discuss them in the 
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context of the appropriate case. In particular, the correspondence will be as 
follows : 


Type trajectory 

Case 

A 

1 (lateral equations) 

B 

1 (normal equations) 

C 

4 


The remainder of the seven possible cases can be constructed from the results 
for these basic trajectory types. 

To clarify the interrelation of the various cases, a flow graph of the 
computations to be discussed in the following sections is shown in figure 1. 
The rectangles contain the variables that are to be computed according to the 
appropriate equations to follow. The indicated tests are made to determine 
whether \p 9 y, and V are required to change. Each possibility, representing 
one of the seven cases, leads to a different path on the flow graph. Dotted 
lines indicate a flow of information in both directions between rectangles in 
a manner made clear by the appropriate equations. The output of the guidance 
generator is indicated by the seven possible outputs and by a list of the 
output variables. 

Case 1. Lateral and Normal Accelerations 

In this section the development of the equations necessary to generate 
the required outputs from the available inputs will be considered. The gen- 
eral case, in which the aircraft is turning and the inclination angle is 
changing, will be considered here, and special cases will be dealt with later. 
Reference should be made to figure 1 for the flow graph clarifying the compu- 
tations to be discussed. 

Lateral equations (Type A trajectory : arc with free end points)- Let 

the points which are selected on the graphic map display of the horizontal 
plane be numbered by Xq , X^ , • • » X^ • • • Xq where 



and the components and y^ represent the coordinates of the ground- track 

display. Figure 2 illustrates the general situation of the three-dimensional 
trajectory which is based on the three chosen points X a , Xg, and Xy . The 
ith and kth segments are straight-line segments and the j th segment is a 
circular segment which forms the smooth transition between the ith and kth 
segments . 

Now consider how the quantities m^ and m^, required for the footprint 
program, can be determined from the three general points shown in figure 2. 

In general, the two points ^ and Xg, connected by the ith line, will be 
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separated by a distance d-^ and the line connecting these two points will 
have a direction cosine SL± given by 


d i 2 ■ «B - V T(X U - V 


( 2 ) 




i 



(3) 


Furthermore, any point X on the ith line connecting X a and Xg is given 
by 


X = 


X a + 


si. 


(4) 


where s is the distance between X a and X. In particular, for the tangency 
points Xp and Xq shown in figure 2, we have 

Xp = X a + s i £ i = X 3 ' (d i " s i )£ i (5 


X q = % + c e £ k 


( 6 ) 


The center of the circular arc X c is given in terms of X p and Xq by a 
relation similar to equation (4) : 


x c = Xp + RAp 


(7) 


X 


C 


X q + R*q 


( 8 ) 


where R is the radius of the circular arc and £ p and £q are the direction 
cosines of the radii between X p and X c » and Xq and X c , respectively. The 
determination of these three quantities is deferred temporarily. Now, by 
equating (7) and (8) and substituting for X p and Xq from equations (5) 
and (6) we obtain 


A 



D 


(9) 


where 


A = (l ± -fc k ) 

D = d.A. + R(A q - A p ) 

The solution is 



( 10 ) 


( 11 ) 
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and these are the quantities necessary to determine the segments required for 
the footprint program. That is, the ith and kth straight-line segments 
shown in figure 2 are 


”1 = s i - c a 

(12) 

m k = s k - c 6 

(13) 


The sequence of calculations in equations (11)- (13) needs some clarifica- 
tion. At the start, when the aircraft is taking off, c a - 0, so that 
from equation (11) is all that is needed in equation (12) to determine m^. 

For the next segment, m^ in equation (13), the s^ comes from a second 
application of equation (11) using three new points starting with Xg while 
the eg is known from the first application of equation (11) . In similar 
fashion the process continues. On the last segment, we must have s^ = dk so 
that 


= d^ - Cg (last segment) (14) 

As for the jth segment, both the subtended angle 0j and the radius R 
are required for the footprint program and for the preceding equation (11) . 

The subtended angle is given by the fundamental relation 

cos 0 j = £ p T £ q (15) 

The footprint program does not recognize a sign on this angle so that the 
input to the footprint program is 

0j = I cos -1 £ph q l (16) 

The sign information is accounted for by specifying R as >0 or <0 
depending on whether 0j was clockwise or counterclockwise. The manner in 
which this is done will be discussed shortly. 

The magnitude of the radius must be determined by the acceleration to be 
allowed. The component of the velocity in the ground- track horizontal plane 
for the jth segment is Vj cos Yj. Hence, the lateral acceleration a^ is 
given by 

(V-j cos Y.)2 

a* = Vj* cos Yj = R (17) 

R should be based on the total acceleration, rather than determined from this 
equation; for this reason the determination of R is deferred temporarily. 

It still remains to determine the equations for £p and &q . It would 
appear that because of the perpendicularity of the radius and the ith and kth 
segments, the conditions 
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( 18 ) 


h\ = 0 

*k\ = 0 ( 19 > 

would be sufficient to define £p and &q. However, this is not the case. In 
the approach to be used here, use will be made of sketch (c) . 

v I 



Sketch (c) 


We will determine whether the extension of line i lies to the left or right 
of line k. Let the line i be extended beyond point Xg . Any point on this 
extension a distance s from X a is given by 


X = X a + sl ± 


or 



Letting x = Xy, we have 



The only unknowns here are y and s and the solution for y is 


y 


y + 

; a 


*y 


- X^ 




( 20 ) 


( 21 ) 


Now it is clear from the geometry of the situation, recalling from equa- 
tion (3) that 5 cos 6^ (see sketch (c)), that 


> 0 and 


yy > y =* left turn 
y Y < y =► right turn 
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< 0 and 


y Y < y 
y Y > y 


left turn 
right turn 


A left turn implies that R lies on the left of segment i while 
turn implies that R lies on the right of segment i. Define the 
of £ p and £ q by 



Now for a left turn we have 



In a similar way, the components of £ q are related to the segment 
we have 



left turn =* < 



For a right turn, the relations for £ p are 



The components of £ q are related in similar fashion to the segment 



right turn =* < 



a right 
components 


k • Thus , 


k . Thus 
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The preceding results 


are condensed by 


the following: 


£ 


P 


£ 


q 



sgnfo^Cy 


sgn[a ± (y 



The orthogonality requirements of equations (18) and (19) are evident. 


( 22 ) 


(23) 


Normal equations (Type B trajectory : arc with constrained end points )- 

Consider now the motion and normal acceleration required to change the incli- 
nation angle. For this purpose the new diagram shown in figure 3 will be 
convenient. Here the attitude is depicted as a function of the distance a 
traversed in the ground track, and the point corresponding to in 

figure 2 is defined by 


z » 5 0 

to reflect the correspondence. Three such points are shown in figure 3. If 
the normal acceleration occurs between the same points Xp and Xq as for the 
lateral acceleration, it is clear that the transition trajectory due to normal 
acceleration is one with fixed end points. Thus, this situation is an example 
of a type B trajectory with constrained end points. 


First, let us determine the transition trajectory between the segments. 
The ith segment has the direction cosine given by 



(25) 


where the value is known, that is, specified; similarly for Ifc . It is 

important to note here that in the equations to follow, y is taken to be 
positive for angles above the horizontal and negative below. Thus, the range 
of interest is -tt/2 < y < tt/2. The point Z p is given by inspection of 
figure 3 as 


Z 


P 


z 


q 



(26) 


(27) 
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where Az is the increase in z along segment j . The curvature required in 
the transition can be found by locating the center Z c of the circular arc as 
follows: 


z c 

Z c 


+ v 

(28) 

+ £qP 

(29) 


where l p and Z q are the direction cosines of the radii shown in figure 3, 
and are represented by 


Z 


P 


Z 


q 



(30) 


(31) 


The determination of the direction cosines will be deferred temporarily. 
Equating equations (28) and (29) and utilizing equations (26) and (27) we 
obtain 


/°a + 

\ z + 

\ a 


s i 



The solution for the unknowns p and Az is 


P 


R0 j 


C 


q 


(32) 


n p ‘ n q 

Az = p(n 0 - n q ) = Re.? ~ — . • (33) 

5 p c q 

and these will be used in the following. It is of interest to note that in 
the limit as ** Yj[, Az = R0j tan Y^, and by inspection of figure 3, this 
can be seen to be the* correct limiting value. 


The radius p just determined can now be used to determine the normal 
acceleration on the j th segment: 


a 


Y 


V V (? P - g q } 

P - ROj 


(34) 
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where equation (32) has been utilized to eliminate p. This equation will be 
used in the next section to determine the lateral turn radius. 


The value of Az determined by equation (33) is used to establish a 
value for the inclination angle during the transition between straight-line 
segments. This matter warrants further discussion. When the inclination 
angles on consecutive straight-line segments are different, that is, Yj_ / Y^, 
the inclination angle Y j on the circular connecting segment will vary in a 
continuous manner between Y ± and Y^. However, the noise footprint program 
only accepts constant Y values. Consequently, in order to meet this 
requirement of the footprint program, it is necessary to determine an appro- 
priate constant value for Yj. This should be done so that the end point Z 
is the same value that would be obtained from the circular segment. Thus, ^ 


Az 

U = ^eT 


n P - n q 
?p - C q 


where equation (33) has been used to obtain the expression on the right side. 

It is still necessary to determine the direction cosines Zp and Zq , 
which have been tacitly assumed to be known in the preceeding development. 
Since p in figure 3 is perpendicular to both the ith and kth segment, we 
have 

= 0 (36 

l k T l q = 0 (37 


Using equation (25) in the above equations yields 


cm:;; 


3-0 


The orthogonality relation is not sufficient to determine the signs, and more 
information is needed. 


Consider now the determination of the signs of the direction cosines Zp 
and Zq. From geometric considerations in figure 3 and the way in which Y 
is defined, it is clear that 

Y i > Y^. =* p lies below the path 

Y^ < Y^. =* p lies above the path 
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as the aircraft proceeds along the trajectory. Thus for Y^ > Y^, 

Cp = cos ^ Y i " l) = sln Y i = n i 

Up = sin^ - f ^ = -cos Y ± = -Z ± 

In similar fashion. 



For the other case in which < Y k , the components of Z p are 



and the components of l ^ are 

? q = _T1 k 
= C k 


The equations in this paragraph are valid for -ir/2 < Y < ir/2. The range of 
Y is restricted because for large angles near 2ir the above sequence of 
equations breaks down, and because [-tt/2, tt/2] is the range of practical 
interest . 


The foregoing can be summarized by the following 


l 


P 



[sgn (y ± 



l 


q 



[sgnCYi 



(40) 


(41) 


where -tt/2 < Y < tt/2. The orthogonality requirements of equations (36) 
and (37) are clearly satisfied. 


Lateral turn radius- The determination of the lateral turn radius R has 
been deferred until now so that it could be based on the total acceleration. 
The total acceleration a Q is given by combining equations (17) and (34): 
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( 42 ) 


a 


2 

o 



+ a Y 


2 



where the j subscript has been suppressed because the velocity is constant 
on segments i, j, and k. Equation (42) can be solved for R for a specified 
value of a Q , and the result is 


R 2 



(V cos Y j ) 4 



(43) 


Since Yj varies from Y j_ to Y^ corresponding to the ends of the j th seg- 
ment, the maximum acceleration will occur at either end, and a maximum opera- 
tion is required as shown in equation (45) below. Further, the footprint 
program requires a sign convention on R to indicate clockwise or counter- 
clockwise turning. The correct sign is obtained from the following: 


sgn(a i ) (y - y y ) > 0 =* 
sgn(a ± ) (y - y y ) < 0 =► 
as can be seen from equation (22) and 

R = sgn[(a ± )(y - y Y )] max — 

n=i ,k a o 


clockwise =* R > 0 
counterclockwise =* R < 0 
the discussion leading up to it.^ 

1/2 


( 


(V cos Y n ) 4 + 


- 


’p 


(44) 


Thus , 


(45) 


This result determines the lateral turn radius such that the total accelera- 
tion will not exceed a specified value. This result is also needed to com- 
plete the determination of the normal equations which have been dependent 
on R. 


Case 2. Lateral Acceleration Only 

This is a trivial special case of Case 1, and the equations given there 
are sufficient. In Case 2, there is no change in the inclination angle 
between consecutive segments, and this includes spiral type motion. The flow 
graph of the computations involved here are shown in figure 1. 


Case 3. Normal Acceleration Only 

The guidance equations for the case in which only normal acceleration is 
required on the transition between straight-line segments are considered here. 
This includes the planar so-called " two- segment" paths. The trajectory 
involved is Type A, that is, an arc with free end points. For this reason. 
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the results for this case are not obtainable from the normal equations for 
Case 1 since the normal equations there were for a Type-B trajectory, an arc 
with fixed end points. That is, the junction points of the transition segment 
for Case 1 were determined by the lateral motion and the normal motion was 
"fit" to these junction points. In the present case the lateral motion is 
absent so that the junction points for the normal equations are free. Since 
the present case is a Type-A trajectory, it is basically the same as for the 
lateral equations of Case 1, and the results there are immediately applicable. 
However, there are minor differences, and in order to avoid confusion, the 
equations and diagram of the situation will be given. Reference should be 
made to figure 1 for a flow graph of the computations . 


Consider the generation of the segment lengths required^by the noise 
footprint program. For the case being considered, we have i/j = 0 and Y ^ 0 
at the junction between straight-line segments. This situation is depicted in 
figure 4. Since the equations for this case will be similar to those devel- 
oped earlier for the lateral case, the distance variables have been italicized 
to distinguish them from those shown in figure 2. The junction points Zp 
and Zq are given by 


z p = 


3 ih = 


- (d. i - s i )Z i 


(46) 


Z q Z 3 + C $\ 


(47) 


where s ^ and C'q are the quantities to be determined. The center of the cir- 
cular arc Z c is given in terms of Zp and Zq by 

Z c = Z p + pZ p (48) 


^c ^q + P^q 


(49) 


where p is the radius of the circular arc, and Zp and Zq are the direction 
cosines of the two radii shown in figure 4. Equating equations (48) and (49) 
and utilizing equations (46) and (47) we obtain 



where 

B = (7. ~ly) 


(50) 


E = d i l i + p(Z q - 7 p ) 


Of the various quantities involved here, Zj^ and ly are still given by equa- 
tion (25), Zp and Zq are given by equations (40) and (41), and d$_ is given 
by 
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(51) 


d l 


d i 


cos Y-£ 


d i 


where d-^ is defined in equation (2). The remaining variable, the radius p, 
must be based solely on the normal equations since there is now no lateral 
motion for this case. For a given choice of normal acceleration, ay , p on 
the j th segment is given by 


P 



(52) 


The solution of equation (50) is 



(53) 


and these are the quantities that are needed in order to determine the segment 
lengths required by the footprint program. The lengths of the ith and kth 
straight-line segments are then given by 


m± = 


o 


a 


(54) 


m k = s k - *3 


(55) 


and they are determined sequentially in conjunction with equation (53) much as 
has been described earlier in connection with equations (12) and (13) . Since 
the footprint program requires segments in the horizontal plane, the ith and 
kth segments to be used in the footprint program are given by 

m ± = m ± cos y ± = rn ± K ± (56) 

m k = m k cos Y k = ^k^k (57) 

The equations simplify on the last segment where the Y*s on consecutive seg- 
ments are the same. Thus if the segment after the kth segment has the same 
inclination angle, it is clear from figure 2 that equation (53) is simply 


Then equation (55) becomes 
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and equation (57) becomes 


m k ~ Wk ~~ c g) cos ^k 


= d k " c 6 cos Y k (58) 

Here the is known from equation (2), and the eg has already been deter- 

mined. The length of the jth segment is, by observation of figure 4, 

nij = Cg(cos Y i + cos Y k ) 

= + 4 k ) ( 59 ) 


This completes the determination of the segment lengths mi, m j , and t% 
required by the noise footprint program. 


The inclination angles Yi, Yj , and Y^ for the ith, jth, and kth seg- 
ments are also required by the noise footprint program. The values for Yi 
and Y^ are known because they were chosen, and the value of Yj varies con- 
tinuously between these values. Since a constant value is required by the 
noise footprint program, a suitable value which gives the correct end point 
Z q is 


fan 


Cg(sin Y^ + sin Y k ) 


m A 


n i + n k 
*i + ? k 


(60) 


This equation can be shown to reduce to the result 


y j - 


Y i + Y k 


(61) 


The remaining variable needed from the guidance generator is the radius 
R which for this case with ip = 0 is 

R = 0 

This condition can be sensed as indicated in figure 1 by determining if con- 
secutive direction cosines, and are the same (or within some small 

value e). Thus, 

|&i - £ k | < e =* R = 0 (62) 


Another variable, useful for later purposes, is the angle 
given by 


8 which is 


cos 


e = i l 

p q 


(63) 
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From equation (63) it can be shown, as expected, that 

e = Y k - y ± (64) 

Case 4. Longitudinal Acceleration Only (Type C Trajectory: 
Longitudinal Acceleration with Fixed End Points) 

Case 4 typifies the Type-C trajectory in which the path is a straight 
line in space with longitudinal acceleration between fixed end points. The 
fixed end points are a result of specifying the points on a land-use map of 
the terminal area between which the longitudinal acceleration is to occur. A 
sketch of the situation is shown in figure 5, where dj is the segment on 
which the acceleration occurs. A flow graph of the computations is shown in 
figure 1. 

Of the variables to be provided to the footprint program, the segment 
lengths, d^, d j , and d^ are given by equation (2). The variable Yj is con- 
stant at the chosen value, and since there is no lateral motion, R = 0. This 
latter condition can be sensed by equation (62) . 


The acceleration between the fixed end points needs to be determined to 
ensure its reasonableness and to enable curved transition trajectories to be 
determined for other cases. For purposes here, we will estimate the accelera- 
tion by determining the constant value of acceleration required to meet the 
boundary conditions, that is, the initial and final velocities, and the dis- 
tance traversed. Analysis shows that a point mass moving in a straight line 
between two points separated by a distance A starting with velocity and 

ending with velocity requires a constant value of acceleration given by 


a v = 



(65) 


To apply this to the present case, we need only to note from figure 5 that 
A = dj/cos Y. Thus, 


(v k 2 - V ± 2 )cos Y 

“v 


( 66 ) 


Case 5. Lateral and Longitudinal Accelerations 

In Case 5, the lateral motion is considered first and the longitudinal 
motion is "fit" within the allowable space. Thus, it is a combination of a 
Type-A trajectory for the lateral motion and a Type-C trajectory for the lon- 
gitudinal motion. The lateral equations are generally the same as for the 
lateral equations developed in Case 1. However, the equation for determining 
the turn radius will be modified as shown below. The longitudinal equations 
are the same as developed in Case 4. A flow graph clarifying the computations 
is shown in figure 1. 
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The lateral turn radius is based on the total acceleration consisting of 
the lateral and longitudinal components and is given by 


„2_ a 2 + _2 

a o a \fj + a V 


pCVj cos Y) r 

2 

_i_ 

’(Vk 2 - v ± 2 ) 

L * J 

T 

2A 


(67) 

where we have used equations (17) and (65) . In this case the distance A can 
be expressed in terms of the lateral turn radius R by 



cos Y 


( 68 ) 


Using the above relationship, we can then solve equation (67) for the lateral 
turn radius. In addition, since the velocity varies from Vf to Vj , it will 
be desirable to ensure that the desired value of total acceleration a Q is 
not exceeded. Further, the sign of R is given by equation (44). Combining 
these ideas, one obtains 


R = 


sgn[ (a 1 ) (y - y ) ] max 
n=i,k 


_1_ 


(V n cos Y) 4 + 


[ 


(v k 2 - v i 2 )cos y 


26 j 


rf 


(69) 


Case 6. Normal and Longitudinal Accelerations 

Case 6 is similar to Case 5 except that the lateral motion is replaced by 
the normal motion. Thus this case is also a combination of a Type-A trajec- 
tory for the normal motion and a Type-C trajectory for the longitudinal motion. 
The normal equations are obtained from Case 3, but the equations for the 
radius must be modified as shown below. The longitudinal equations are the 
same as developed in Case 4. A flow graph of the computations required for 
this case are summarized in figure 1. 


The radius p in this case is determined from the total acceleration 


a 2 =a 2 +a 2 
a o a Y a V 


/V \ 2 / v k 2 - V \ 2 

{—) — ) 


where we have used equations (52) and (65) . Now the distance A is 

A = p0 


(70) 


(71) 
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Combining the latter two equations and solving for 


p, we obtain 


P 


= max 


n=i,k a o 



( 72 ) 


Case 7. Lateral, Normal, and Longitudinal Acceleration 


Case 7 combines the motion due to accelerations in all three axes- Thus 
it is a combination of Types A, B, and C trajectories. The lateral and normal 
equations are generally the same as for Case 1. However, the turn radius will 
be modified as shown below to account for all the components of acceleration. 
The longitudinal equations come from Case 4 as indicated below. The flow 
graph of the computations for this case is shown in figure 1. 

The turn radius R is determined from the total acceleration which is 


a 2 = a 2 
a o a ip 


+ a,/ + a 


V 


(73) 


The first component is given in terms of R by equation (17) and the second 
component is given by equation (34). The third component is given by equa- 
tion (65). However, it is difficult to determine the distance A in that 
equation because of the combined motions. We will be conservative and take 
the distance traversed due to lateral acceleration. Thus the distance A is 
given by equation (68) where the Y is now Yj from equation (35). That is, 


A 



cos Yj 


(74) 


Hence we have 



-.2 

COS Yj ) 2 

~R 



i 


(V k “ V i ) 2 cos Y 


2R0 


I 


To suitably limit a Q it is necessary to know where the maximum occurs on the 
jth segment. If Vj and cos Yj are greatest at the same end of the jth 
segment, the maximum acceleration will occur at that point. If Vj and cos Yj 
are not maximum at the same end of the jth segment, the maximum acceleration 
may theoretically occur anywhere on the segment, and it is difficult to deter- 
mine where. However, the maximum will nearly always occur where Vj is 
largest since changes in cos Yj are inherently very small. Thus we will 
base the turning on the acceleration at either end of the jth segment. 
Further, the sign convention of equation (44) will be required. Thus we have 
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( 75 ) 


R = 


sgn[ (a. ) (y - y y )] max 
n=i,k 


1 



I 


(V n cos Y n )‘ 


2 i/2 

- Videos Y n l 1 

29 J JJ 


AIRCRAFT CONTROL GENERATOR 


The generation of the optimal aircraft controls is the subject of this 
section. Given the guidance information, it is desired to have the aircraft 
follow the commanded trajectory by making the proper choice of the available 
controls. One of these controls, the engine parameter information, will be 
required by the noise footprint program. The solution will depend on the 
class of aircraft to be considered, that is, whether the thrust vector is 
rotatable and whether the thrust affects the aerodynamic forces. 

The aircraft equations of motion will be needed as a starting point for 
determining the controls. An examination of the equations of motion in 
various coordinate systems shows the velocity coordinate system to be the most 
appropriate. A set of orthogonal coordinates is chosen to be in the direction 
of the velocity vector, perpendicular to the velocity vector in a vertical 
plane, and perpendicular to the velocity vector in a horizontal plane. A 
diagram which displays the various forces involved is shown in figure 6. The 


point-mass dynamic equations are 

mV - T cos (a + q) - mg sin Y - D (76) 

mVY = [L + T sin(a + n)]cos <j> - mg cos Y (77) 

mV^ cos y = [L + T sin(a 4- q)]sin <j> (78) 

The kinematic equations are 

z = V sin Y (79) 

x = V cos Y cos ip (80) 

y = V cos Y sin ^ (81) 

The aerodynamic forces are given in general by 

D(V,a,T,6 F ) = pV 2 SC D (a,T,6 F ) (82) 

L(V,a,T,6 F ) = -| pV 2 SC L (a,T,6 F ) (83) 
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where the inclusion of T represents the effect of thrust on the aerodynamic 
forces. These equations of motion have three state variables, V, Y, and 
and four control variables, T, n, <j>, and a which are to be determined. It is 
worth noting that the equations in this coordinate system seem to represent a 
good balance in simplicity between the dynamic and kinematic equations. For 
example, if the dynamic equations are written in an earth- fixed coordinate 
system, the kinematic equations become simpler but the dynamic equations are 
made more complicated, and a solution for the optimal controls is difficult to 
obtain. 


Now consider the solution of the equations of motion for the controls. 
Equations (76) to (78) become, when slightly rearranged. 


T cos (a + n) = mg sin Y 4- D + mV 

T . m . / , x mg cos Y + mVY 

L + T sin(a + n) = — ® 7 

cos <p 

T , m . / , cos Y m(V cos Y ) 2 

L + T sm(a + n) = 7— L = — ^ — 

sm <p R sxn $ 


(84) 

(85) 

( 86 ) 


where in the latter the following basic relationship for the lateral accelera- 
tion a^ has been used: 


The control for 
and ( 86 ): 


a^ = Vip cos y 


(V cos Y ) 2 
R 


(87) 


can be obtained uniquely by equating equations (85) 


= tan' 


-1 cos Y 

g cos y + Vy 


= tan 


-1 


(V cos Y ) 2 
R(g cos y + Vy) 


( 88 ) 


The controls T and a can be related by solving equations (84) and (85) , and 
the result is 

T 2 = (mg sin Y + D + mV ) 2 + ( ™ g COS Y + mV ^ - l\ (89) 

1 cos 9 I 

The two controls n and a can be related by means of equation (84) 


n 


cos 


■( 


- if mg sin Y + D •+ mV 


9 - 


(90) 


These three equations comprise the relationships to be satisfied for the four 
controls . 


It is now necessary to examine the conditions under which T is constant 
since the footprint program will accept only constant values. Sufficient con- 
ditions are 
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Y = V = 0 
and ifj = constant ) 
Thus, the control equations are 


( 91 ) 


= tan 


-1 Vjfe. = tan -l V... C Q 8 Y 


Rg 


r 2 = [mg sin Y + D(V,a,T,5 F )] 2 + - L(V,a,T,6p)J 

• l/mg sin Y + d\ 

r - * — ) ■ ° 


n = cos 


(92) 

(93) 

(94) 


If the engine parameter T is not constant, it is approximated here by a 
linear variation between the two constant values of the preceding and follow- 
ing segments. The footprint program contains this capability. Let and T k 
be the constant value of the engine parameter on the ith and kth segments 
and let Tj be the non-constant value of the engine parameter on the jth 
segment. Then if one sets Tj = T k the footprint program will cause Tj to 
vary with the distance a as 


Tj (a) = T ± + (T k - T i )0(a) 


where 


6(a) = 3 — _ 0 < a < 304.8 

This function causes Tj to vary linearly from Tji to T k in a fixed dis- 
tance of 304.8 m (1000 ft). A better procedure would be to vary Tj over the 
actual computed distance between the ith and kth points. This could be done 
from equations given earlier, but would require modification to the footprint 
program. 

If the conditions for constant engine parameter T are met, it will be 
necessary to solve equations (92) to (94) . The solution of these equations 
depends on the class of aircraft to be considered and there are several possi- 
bilities. The thrust vector angle rj may or may not be rotatable and the 
thrust vector may or may not affect the aerodynamic forces involved. Thus in 
the following sections, we will discuss three classes of aircraft. 


Class A Aircraft: Thrust Angle Fixed 

The Class A aircraft is typified by the CTOL aircraft in which the con- 
trol n, the angle of thrust vector, is fixed with respect to the aircraft. 
For simplicity, n is taken to be zero. An aircraft of this class typical of 
the current fleet is described in reference 3. Further, the aerodynamic 
forces, lift, and drag are not affected by the engine thrust so that they are 
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characterized by D(V,a) and L(V,a) for a fixed flap setting. Thus the con- 
trol equations (92) to (94) become 


. -i Vip . -i cos Y 
+ - tan — - tan 

- [mg sin Y + D(V,a)] 2 + Y - L(V,ol)1 


T = mg sin Y + D(V,oQ 


cos a 


(95) 

(96) 

(97) 


The first equation gives <j> based on information from the guidance generator. 
The second and third each express a relationship between the two controls T 
and a, so that there is a unique solution. We will not deal with various 
esoteric approaches for solving these equations. Here we will simply solve 
the latter two equations simultaneously by a one-dimensional search over a 
such that equations (96) and (97) are equal. This procedure will be illus- 
trated by a later example. 


Class B Aircraft: Thrust Angle Variable 

For the Class B aircraft the angle of the thrust vector is rotatable, 
thus providing another degree of control. In this case equations (92) to (94) 
form a redundant set of the form 



cj) = f^(states) 



(98) 


T 2 = frp(a,T,(j), states) 



(99) 


D = f^ (a, T, states) 



(100) 

These equations contain four control functions, <j>, T, a, and q 
states, V, y, and ip . The redundancy is clear because T and n 
determined from equations (99) and (100) for any arbitrary a. 

and three 
are uniquely 
To remove this 


redundancy requires an additional condition. We will require the additional 
condition that the noise be minimized by minimizing the engine thrust T. For 
this reason, interest centers on either equation (93) or (99) in order to 
satisfy the additional condition. The other control equations are not 
involved with this additional condition because the determination of the con- 
trol <J> from equation (92) or (98) is . independent of the other control 
variables, and the determination of r\ from equation (94) or (100) follows 
directly from the solution for the controls T and a. Thus, the additional 
condition to be imposed is 

Min T: T 2 = f T (a ,T, <J>, states) (101) 

a 


subject, of course, to any particular constraints on a and n» depending on 
the type of vehicle. How this equation is solved depends on the aerodynamic 
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characteristics of the aircraft, that is, whether the aerodynamic forces are 
affected by the engine parameter, either unavoidably or purposely. Thus, 
there are two principal types to be considered and two distinctly different 
approaches for a solution. 

Class B1 aircraft: aerodynamics independent of engine parameter - For 

certain classes of V/STOL aircraft, the aerodynamic forces, lift, and drag, 
are not affected by the engine parameter T. These Class B1 aircraft have 
been widely used in various studies such as that of reference 4. Some lift- 
fan aircraft are in this category and many design studies have been made for 
these aircraft (refs. 5 and 6). The lift and drag for such vehicles are 
functions of only two variables; that is, D(V,a) and L(V,a) . The consequence 
of this is that the minimization in equation (101) becomes of the simpler form 

Min T: T 2 = f T (ot , <f> , states) (102) 

a 

If the aerodynamic forces are given in tabular form, the minimization in 
equation (102) can be solved by a simple one-dimensional numerical search 
over a . 

If the aerodynamic forces can be approximated by a simple expression, the 
minimization can be accomplished analytically. The well-known Kuhn- Tucker 
problem is to 


Max f(a) 


subject to 


a < a L 


where is the maximum permitted value of a. Since we want to minimize T, 

we set f(a) = -T 2 where T 2 is given by equation (93). The approximations 
to be used for lift and drag are 


D(V,a) = j pSV 2 C D = j pSV 2 (D 0 + DjCi 2 ) 

L(V,a) = \ pSV 2 C L = \ pSV 2 (L 0 + I^a) 

Using these approximations, we obtain 

[ 1 -| 2 1 
mg sin y + pSV 2 (D Q + DjOt 2 ) 

= -(a + ba 2 ) 2 - (c - da) 2 


(103) 

(104) 


(105) 
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where a, b, c, and d are functions of Y, V, and (p ; the functions are 
apparent from the above equation. The inequality constraint can be made into 


an equality constraint by the slack variable y: 

a + y 2 - a L = 0 (106) 

Thus, the Lagrangian for this problem is 

L = -(a + ba 2 ) 2 - (c - da) 2 - A(a + y 2 - a-jO (107) 

and the necessary conditions for a minimum are 

L a = -4b 2 a 3 - 2 (2ab + d 2 )a + 2cd - A = 0 (108) 

L y = -2Ay = 0 (109) 

= -(a + y 2 - a L ) = 0 (110) 


From equation (109) it is seen that there are two possible candidates for 
the optimum. Either A = 0 or y = 0, and it is necessary to examine both 
possibilities to see which will lead to the solution. If A = 0, then a is 
not on the boundary and the optimum a is obtained from the solution of the 
cubic equation in a, equation (108) with A = 0. On the other hand, if 
y = 0, the optimum a is on the boundary. In either case the solution for 
a can be used to determine T from equation (105) or (93) and the thrust 
angle r\ by equation (94). A later example will illustrate this discussion. 

Class B2 aircraft: aerodynamics dependent on engine parameter - For many 

V/STOL aircraft the effect of the thrust parameter on the aerodynamic forces 
cannot be neglected. Thus the lift and drag are functions of three variables: 
L(V,a,T) and D(V,a,T). There are many examples of Class B2 aircraft, such as 
the externally blown flap, the augmentor wing, and the upper surface blowing 
jet flap aircraft (refs. 7-10). 

For this class of aircraft, the thrusting force is a function of a con- 
trol variable 0, such as throttle setting 

T = f T (0) 

and the aerodynamic lift and drag coefficients are functions of two variables 
for a fixed flap setting: 


C L = f L (o,0) 

Cp = f d 9 ® ) 

Generally, Cp and Cp will be expressed in graphic form. The variable 0 in 
f L and fp reflects the interaction of thrust with the aerodynamics, whether 
it is unavoidably introduced or whether it is purposely introduced. The con- 
sequence of this interaction is that the minimization in equation (101) is 
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difficult because T appears on both sides of the equation and because Cl 
and C D contained in frp will generally be in graphic form. 

The minimization in equation (101) can be solved by a two-dimensional 
search for the optimal values of T and a. Although various sophisticated 
approaches for solving equation (101) could be utilized, here we will simply 
do the following. First, define e 2 to be the square of the error by which 
equation (99) or (93) is not satisfied, that is, 

e 2 = T 2 - [mg sin Y + D(V,T,a)] 2 - ~ - L(V,T,a)J (111) 

Second, examine the mapping Txa e by flooding the Txa space. Third, 
select the set (Txa: e = 0}, that is, the set satisfying the equations of 
motion. And fourth, carry out the minimization indicated in equation (101) by 
selecting the member for which T is minimum. More concisely, 

Min{T: T G (T,a) 3 e = 0} 

a 

If desired, the whole process can be carried out directly on the computer so 
that results of intermediate steps are not displayed. 


EXAMPLES 


In this section we will illustrate the results of applying the preceding 
method to the evaluation of noise impact at a typical airport. Figure 7 shows 
an actual airport, including the airport boundary, the runway, and the sur- 
rounding zoning pattern as indicated by the codes. Each of the irregular- 
shaped areas (even with the same zoning) corresponds to a different land 
value, with the higher priced land indicated by planned residential zoning 
generally lying in the southern portion of the map. The noise impact is 
determined by superimposing a noise footprint on the zoning pattern and by 
evaluating the total value of land subjected to unacceptable noise levels. 

For this purpose we will consider the generation of the takeoff footprints for 
the various types of aircraft. The same equations are valid for landings. 

First, we will consider an application to an aircraft of Class A, a 
representative current CTOL aircraft. For this aircraft class the thrust 
angle is fixed and the aerodynamic forces are functions of only velocity and 
angle of attack for a fixed flap position. The data needed here for such an 
aircraft were taken from reference 3 and are as follows: 

W = 781,047 N (175,595 lb) 

S = 144.93 m 2 (1560 ft 2 ) 

C L = 0.60 + 0.1065 a 
C D = 0.0845 + 1. 136xl0 -4 a 2 
T max = 192,154 N (43,200 lb) 

The noise source characteristics, which are needed as stored information in 
the footprint program, consist of data depicting EPNdB as a function of two 
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variables, the slant range to the observer and the magnitude of the engine 
parameter (i.e., thrust level). Representative data for current CTOL aircraft 
are shown in figure 8. 

A facsimile of the scope display of the airport and surrounding zoning 
patterns is illustrated in figure 9. Shown here is the coordinate system 
starting at the origin of the runway. Also shown on this figure is the flight 
path passing through the chosen way-points (indicated by small circles) with 
chosen inclination angles. The way-points are put in directly on the scope 
display by the operator and sensed by the computer, while the inclination 
angles are put in on a keyboard. With this information the guidance generator 
calculations are made in accordance with the equations for Case 3 and the flow 
diagram of figure 1. The computations made by the guidance generator are sent 
directly to the footprint program and are not normally made available to the 
operator. However, for purposes of illustration, these calculations are shown 
in columns 2, 3, and 4 of the following: 


Segment 

Guidance generator output 

Control generator output 

Length, m (ft) 

Y, deg 

R, m (ft) 

T, % 

<f>. deg 

a, deg 

1 

1,524.0 (5,000) 

0 

0 

100 

0 

0 

2 

954.87 (3,132.8) 

3.75 

0 

Varies 
Set to 82.6 

0 

— 

3 

30,004.61 (98,440.3) 

7.50 

0 

82.6 

0 

5.55 


The aircraft control generator calculations are determined by the solu- 
tion of equations (95) through (97). The solution of equations (96) and (97) 
involves a one-dimensional search for the value of a which satisfies both 
equations. This is illustrated in figure 10 for the third segment, the only 
one needed for this simple trajectory. Here the variation of equations (96) 
and (97) with a is shown, and it can be seen that there are two solutions 
close together. The desired solution is the one indicated because it requires 
less thrust and hence results in less noise. The controls for the three seg- 
ments are summarized in the last three columns of the above table. 

The information sent to the footprint program consists of columns 2, 3, 

4, and 5 of the above table. The output of the footprint program is given by 
the single-event noise contour in figure 9. It is seen that the noise impact 
on the surrounding community is substantial. Moreover, the limited maneuver- 
ability of this aircraft relative to the later V/STOL aircraft means that 
little can be done to alleviate this effect. For example, the ground- track 
corresponding to a 0.1-g turn starting at an altitude of 304.8 m (1000 ft) is 
shown by the dotted line in figure 9; the corresponding contour is not shown 
but it is clear that the noise impact would still be substantial because of 
the limited maneuverability of this aircraft. 

The second example is for Class B1 aircraft which are typified by rotat- 
able thrust vectors and aerodynamic forces that are independent of thrust. A 
vehicle that fits this category approximately is the lift-fan 100-passenger 
commercial short-haul transport studied in references 5 and 6. To avoid 
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repetition, we will illustrate for this example only the aircraft control 
equations. For this purpose we will take the flight path of the aircraft and 
the guidance equations to be the same as for the next example, since the air- 
craft for both examples are highly maneuverable and can readily follow the 
required path. This path (shown in fig. 11) has been chosen as a potential 
candidate for minimizing the noise effect on takeoff because it largely avoids 
much of the residential property. It might be noted that this path does not 
require the present lift-fan aircraft to utilize its full capability since the 
path is a STOL flight path whereas the lift-fan aircraft has vertical flight 
capability . 

The aircraft control generator utilizes the output of the guidance gener- 
ator (see results given in the next example) and the basic aircraft data to 
determine the optimal control functions. The aircraft characteristics needed 
for the simplified model here are taken from references 5 and 6 and are 

W = 561,782 N (126,300 lb) 

S = 73.21 m 2 (788 ft 2 ) 

C L = 0.94 + 0.1017 a 
C D = 0.18 + 0.001342 a 2 
T max = 646,050 N (145,245 lb) 

«L = 10 ° 


where the expressions for and Cj) have been fit to aerodynamic data in the 

above references. The permissible magnitude of the angle of attack ctL is 
influenced by several considerations. First, vertical gust considerations 
require that an a safety margin be provided to insure that the peak of the 
lift curve not be exceeded. Reference 11 indicates the magnitude of this 
safety margin should be sin~ 1 (20/V) ° . Second, the angle of attack should be 
limited for reasons of passenger comfort, that is, to restrict the pitch 
attitude of the aircraft. A reasonable limit on a to satisfy both consider- 
ations will be taken to be = 10°. The solution of the control equations 

is required for each segment of the flight and follows previous discussion. 

For segment 3, for example, the control <f> = 0. The determination of the con- 
trol a involves satisfying the necessary conditions given by equations (108) 
through (110) . For this purpose the Lagrangian is 

L = -(23,202.92 + 17.576 a 2 ) 2 - (112,256.93 - 1,331.939 a) 2 - A (a + y 2 - a L ) 

Now if A = 0, Ly = 0 is satisfied. Then 

L a = 1235.64 a 3 + 5.17938xl0 6 a - 2.99039xl0 8 + A = 0 

is satisfied by a = 41.13°. But then = - (a + y 2 - 10) = 0 cannot be 
satisfied. On the other hand, if y = 0, Ly = 0 is satisfied, and = 0 

is satisfied by a on the boundary, that is, a = aL = 10°. Finally, L a = 0 
is satisfied for A - 2.460><10 8 . Thus, all three necessary conditions for the 
optimum are satisfied with a on the boundary. The other controls T and n 
are then found using this value of a and equations (93) and (94). The 
results for all the controls are summarized for all segments in the following: 
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CONTROL GENERATOR OUTPUT 


Segment 

<j), deg 

a, deg 

T, % 

n, deg 

1 

0 

0 

100 

— 

2 

— 

— 

Varies 
Set to 70.3 

— 

3 

0 

10 

70.3 

65.8 

4 

5, 7, 9, 11 

0 

10 

Varies 
Set to 61.7 

61.7 

57.97 

6, 8, 10 

14.2 

10 

64.6 

59.0 


The third application is to a vehicle of class B2 which is characterized 
by a rotatable thrust vector and aerodynamic forces that are significantly 
affected by thrust. In accordance with previous discussion, the characteris- 
tics of such vehicles are represented in sketch (d) where 0 is the primary 
engine parameter controlling thrust. From an input- 
output viewpoint, the internal structure is unim- 
portant. However, in order to generate and interpret 0 ^-j 

the input-output characteristics, it is necessary to 
become concerned with the internal structure for a 

particular aircraft. It will be convenient for this — — 

purpose to use a research version of the augmentor a _^ c 

wing aircraft because of the ready availability of 
the required data. 



L 

D 


The augmentor wing aircraft is provided with Sketch (d) 

direct lift by means of the hot exhaust of the 

engines which can be rotated downward from horizontal, covering a range of q 
from 0 to 98°. In addition, the normal lift of the wings has been augmented 
by the addition of a cold air duct and flap redesign that exhausts the cold 
air along the wing. This results in added lift due to increased circulation. 
For such a vehicle the internal structure is as shown in sketch (e) , where 0 


is the engine throttle angle, T is 
the thrust due to the hot exhaust gas, 
T q is the cold thrust, and Cj is 
the cold thrust coefficient defined by 
Cj = T^/qS. The characteristics of 
the unlabeled boxes were obtained 
from reference 9 and are shown in 
figures 12 and 13. 

It is worth reiterating that the 
aerodynamic characteristics shown are 
representative of this class of air- 
craft in which thrust interacts with 
aerodynamic forces. That is, the 
lift and drag coefficients are 



I J 

Sketch (e) 
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invariably given as a function of two variables, the angle of attack and an 
interaction parameter (such as Cj) related to the thrust (see data given in 
ref. 7 for a variety of configurations). Although the physical mechanism of 
this interaction depends on the specific aircraft type, the input-output char- 
acteristics are all described by similar kinds of data. 

The remaining characteristics of the vehicle under study are (see ref. 9): 

W = 177,920 N (40,000 lb) 

S = 80.36 m 2 (865 ft 2 ) 

T max = 61,605 N (13,850 lb) 

Finally, the noise characteristics may vary depending on the aircraft type and 
the level of noise suppression technology incorporated into the aircraft. 

Since the sensitivity of noise effect to the level of noise suppression is of 
inherent interest, we will take the two different levels, A and B (fig. 14) to 
illustrate the results. Level B is taken from unpublished FAA data which are 
approximately the same as data given in reference 7, while level A has been 
increased by 10 dB in order to investigate the sensitivity to this parameter. 

The flight path chosen for the aircraft to fly is shown by the facsimile 
of the scope display in figure 11. The way-points, which are chosen directly 
on the scope display, are indicated by the small circles, and the inclination 
angles, which are put in through the keyboard, are also indicated. With this 
information, the guidance generator calculations are made in accordance with 
previous discussion. Again, the results of these computations are normally 
sent directly to the footprint program and are not available for inspection. 
However, the results are illustrated in the following. 


GUIDANCE GENERATOR OUTPUT 


Segment 

Y, deg 

R, m (ft) 

V, knots 

Length or angle, 
m (ft) or deg 

i 

0 

0 

Varies 

194.9 (639.4) 

2 

4.75 

0 

70 

218.2 (715.8) 

3 

9.5 

0 

70 

196.5 (644.6) 

4 

9.5 

0 

Varies 

914.4 (3,000) 

5 

9.5 

0 

100 

239.6 (786.2) 

6 

9.5 

+1050.71 (+3447.2) 

100 

38.8° 

7 

9.5 

0 

100 

254.3 (834.2) 

8 

9.-5 

+1050.71 (+3447.2) 

100 

27.4° 

9 

9.5 

0 

100 

847.2 (2,779.4) 

10 

9.5 

+1050.71 (+3447.2) 

100 

45.8° 

11 

9.5 

0 

100 

6037.4 (19,807.8) 


The information in columns 2, 3, and 5 are sent to the footprint program while 
the information in columns 2, 3, and 4 are sent to the aircraft control 
program. 

The generation of the optimal controls follows the theory previously 
described. The mapping Txa e has been investigated numerically and the 
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set obtained for which e = 0 (see eq. (111)). For illustration, this set is 
displayed in figure 15 for segment 5 of the flight path. Also shown is the 
other control q. From this figure, the removal of the redundancy by the min- 
imization in equation (101) so as to minimize the noise amounts to simply 
choosing the minimum T subject to constraints on a and n . It can be seen 
that the true minimum T occurs at essentially the boundary of the q con- 
straint n < 0 as shown. The computer program does not locate this point 
precisely because of the finite increments in T used (limited to about 
80 lb). Thus, the computer program solution is slightly in error. Note that 
this error results in a significant difference in q between the two points. 
However, this is an academic point that is not of great concern because the 
engine parameter T, which determines the noise impact, is essentially the 
same for either solution. A more sophisticated optimization procedure could 
locate the true minimum precisely, but this was not considered necessary here. 
A summary of the results of calculations for all the controls over the entire 
flight path is shown in the following. Of the controls shown, only the data 
for T needs to be sent to the footprint program. 


CONTROL GENERATOR OUTPUT 


Segment 

<f>, deg 

a, deg 

T, % 

n, deg 

1 

0 

0 

100 

— 

2 

3 

0 

10.96 

Varies 
Set to 55.3 

55.3 

10.55 

4 

5, 7, 9, 11 

0 

2.51 

Varies 

Set to 59.78 
59.78 

9.86 

6, 8, 10 

14.24 

2.95 

60.36 

8.5 


The output of the footprint program is displayed in figure 16 by the 
single-event 90 EPNdB contours. The two contours represent the results for 
the two different levels of technology. For the level A technology, it is 
seen that the noise effect on the residential property is fairly small because 
of the capability of the aircraft to fly a curved path that avoids most of the 
noise-sensitive areas. Further, the steep flight-path capability contributes 
to this reduction. On the other hand, the level B technology is sufficient to 
confine the noise effect to within the airport ’boundary . Comparisons such as 
this are of obvious use in investigating the sensitivity of the actual noise 
impact to the noise technology level of the aircraft. 


CONCLUDING REMARKS 


A few points regarding this paper are worth emphasizing. The methodology 
presented here allows complete flexibility in determining the desired flight 
path of the aircraft. This is important because advanced V/STOL aircraft have 
much greater maneuvering capability than the conventional current CTOL 
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aircraft. If the advantages of this greater maneuverability in terms of 
increased airspace capacity, reduced delays, reduced noise, etc., are to be 
shown, the methodology must be flexible enough to be able to account for these 
greater capabilities. The methodology discussed here is flexible enough to 
permit the simultaneous consideration of maneuvers in all three axes. Also, 
the transitions between straight-line segments are uniquely defined, that is, 
there is no necessity to iterate between the transition flight paths and the 
accelerations . 

There are a number of areas for further research and development. First, 
since this report was restricted to the development of the equations, further 
work should deal with the implementation of the computer program of the models 
to provide the desired on-line interactive capability. Second, further devel- 
opment of the footprint program that was used could be made in the following 
directions: 

1. The speed of the program could be increased for better man-machine 
interaction. 

2 . The maximum number of segments allowed could be increased from the 
present four. 

3. A relatively minor change could be made to enhance on-line operation 
so as to give the operator control of the termination of the contour computa- 
tions rather than allowing the program itself to terminate the computations, 
sometimes prematurely. 

4. The program could be changed to avoid the possibility of producing an 
incorrect noise level at a given observer position by basing the computations 
on the segment producing the greatest noise rather than on the closest segment. 

5. The program could be improved to incorporate anisotropic noise 
characteristics . 

Finally, the generation of the aircraft control functions was not dealt with 
adequately. Sophisticated optimization methods could be utilized to great 
advantage to seek out the desired solutions more efficiently than was done 
here . 


Ames Research Center 

National Aeronautics and Space Administration 

Moffett Field, California 94035, November 17, 1977 
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Figure 3.- Illustration for Case 1: Normal trajectory with fixed end points 
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Figure 7.- Typical airport and surrounding community 
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Figure 11.- Flight path for example STOL aircraft. 



THROTTLE ANGLE 0 , deg 

Figure 12.- Hot and cold thrust versus throttle angle for example jet STOL 

aircraft . 
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Figure 15.- Optimal controls for example powered-lift STOL aircraft of 

Class B2; Segment 5. 




Figure 16.- 90 EPNdB noise footprints for example STOL aircraft. 
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